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Abstract 

We construct the two- and three-nucleon potential based on the most general chiral effective 
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a power counting scheme consistent with this projection formalism. In contrast to previous 
results obtained in old-fashioned time-ordered perturbation theory, the method employed 
leads to energy-independent potentials. We discuss in detail the similarities and differences 
to the existing chiral nucleon-nucleon potentials. We also show that to leading order in 
the power counting, the three-nucleon forces vanish lending credit to the result obtained by 
Weinberg using old-fashioned time-ordered perturbation theory. 
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1 Introduction 



One of the most important and most intensively studied problems in nuclear physics is deriving 
the forces among nucleons. At present, we are not able to perform direct calculations of these 
starting from quantum chromodynamics (QCD). Instead of this the alternative approach based 
on the chiral symmetry of QCD was pioneered by Weinberg [|J || and developed by van Kolck et 
al. HHH]. In this method one starts from the most general Lagrangian for pions and nucleons. 
Apart from the usual symmetry constraints, spontaneously broken chiral symmetry is taken into 
account. Consequently, the couplings of pion fields, which play the rule of Nambu-Goldstone 
bosons within this formalism, to nucleon fields are of derivative type. The most general pion- 
nucleon Lagrangian is then given by infinite series of terms with increasing number of derivatives 
and nucleon fields. Explicit breaking of chiral symmetry, leading to the nonvanishing but small 
pion mass, can also be treated systematically. 

The Lagrangian constructed in this way can then be used to derive the expansion of low energy 
S-matrix elements in powers of some low energy scale Q. To be more precise, the 4-momenta 
of the external pions and the 3-momenta of the external nucleons must be no larger but of the 
same order as Q. It is therefore natural to treat the nucleons nonrelativistically. More precisely, 
the dimensionless expansion parameter is Q/A x , with A x the scale of the spontaneous symmetry 
violation. Although the value of A x is not precisely known, it is related to the scale of the higher 
mass states in QCD, i.e. A x can be chosen to be the mass of the rho or 2nF n ~ 1.2 GeV, with 
F n = 186 MeV the weak pion decay constant.^ All other degrees of freedom such as higher 
mass mesons, antinucleons or other baryons are less relevant for low-energy physics and may be 
considered as to be integrated out of the theory. The information about them is then hidden in 
parameters of the Lagrangian. Only the A(1232) resonance, which has a mass close to the one 
of the nucleon and is strongly coupled to the pion-nucleon system, is sometimes kept explicitly. 
Once the effective Lagrangian is constructed, one can formulate the power counting rules in 
order to determine the power of Q for any given Feynmann diagram. It was shown in Q that 
chiral symmetry guarantees a perturbative expansion of low-energy processes because the power 
of Q for any given process is bounded from below. In particular, only a finite number of tree 
and loop diagrams can contribute to given processes at any finite order. 

This scheme works fairly well for purely pionic or processes with only one nucleon 0H§. 
However, complications arise when one tries to treat processes with few nucleons caused by 
infrared divergences (in the static approximation for the nucleons), which destroy the power 
counting. These difficulties are in principle to be expected: the presence of low-energy bound 
states signals the failure of perturbation theory. One way of solving this problem, first introduced 
by Weinberg Q, is to use the time-dependent ("old-fashioned") perturbation theory instead 
of the covariant one. The expansion of the S-matrix, obtained using this formalism, has the 
form of a Lippmann-Schwinger equation with an effective potential, defined as a sum of all 
diagrams without pure nucleonic intermediate states. Such states would lead in "old-fashioned" 
perturbation theory to energy denominators, which are by a factor of Q/mj^ smaller then those 
from the states with pions and which destroy the power counting. The effective potential is free 
from such small energy denominators and can in principle be calculated perturbativly to any 
given precision. Many other equivalent schemes are known leading to same results like that due to 
Bloch and Horowitz || or the Tamm-Dancoff approximation flQ|| . The effective potential, derived 
using the old-fashioned time-ordered perturbation theory, possesses one unpleasant property: it 
is in general explicitly dependent on the energy of incoming nucleons and as a consequence of 
this is not Hermitean. Furthermore, the nucleonic wave functions are not orthonormal in this 
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approach [11]. 



One can avoid these problems by using the method of unitary transformation. It was already 
applied succesfully in cases where one has an expansion in a coupling constant, such as the 



pion-nucleon coupling, see e.g. refs. [12|. Our aim is to show how to apply this method to 
the case of chiral perturbation theory for pions and nucleons, in which the small momenta of 
external particles play the role of the expansion parameter. In fact, our considerations are more 
general since they can be applied to any effective field theory of Goldstone bosons coupled to 
some massive matter fields. 

In section 2 we give a short introduction to the Bloch-Horowitz scheme and the method of 
unitary transformations as pionneered by Okubo [|1|] and by Fukuda et al. flq] . In section 3, 
we show how to apply this formalism to the chiral invariant effective Lagrangian for pions and 
nucleons. We use these ideas in section 4 in order to calculate the two-nucleon force in leading 
and next-to-leading order using the same effective Lagrangian as in Q and compare our results 
with those obtained using "old-fashioned" perturbation theory. We also derive the leading order 
three-nucleon force and shows that it vanishes, a result first found by Weinberg [0. Section 5 
contains the summary and conclusions. Various technicalities are relegated to the appendices. 

2 Bloch-Horowitz scheme and the method of unitary transfor- 
mation 



A good introduction to both schemes is given in [11]. We briefly repeat the main points, mostly 
to establish our notation and keep the paper self-contained. Furthermore, we need the basic 
equations to develop the chiral power counting scheme in the following. 
We start from the Hamiltonian 

H = Hq + Hj (2.1) 

for pions and nucleons with Hq (Hj) denoting the free (interaction) part. The whole problem for 
an arbitrary number of interacting particles can be cast in the form of a Schrodinger equation, 

H\V) = E\&) . (2.2) 

In order to solve this equation for nucleon-nucleon scattering, it is advantageous to project it 
onto a subspace \<j>) of the full Fock space \^>) which contains only nucleonic states. We shall 
denote the remaining part of the Fock space by |\&) = \<p} + \ip). Let r/ and A be projection 
operators on the states \cj>) and \ip) which satisfy rf = 77, A 2 = A, tjX = Xr/ = and A + rj = 1. 
Eq. (|2.2| ) can now be written in the form 

( n H n V HX\( \4>) \ _ ( \<j>) \ , , 

{ XHr, XHX ) { m ) ~ E { W J • (2 ' 3) 

In the Bloch-Horowitz scheme § one obtains from the second line of this equation an expression 
for the state 

H) = E _\ HX XHr lW ( 2 - 4 ) 
which can then be used to reduce the first line of eq. (|2 .3|) to 



(H + V e s(E))\<t>)=E\cf>) (2.5) 
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with an effective potential V e ff(E) given by 



V cS (E) = r)H lV + r lHiX E _ XHX XH lV . (2.6) 



Expanding the denominator in eq. ( |2.6[ ) in powers of Hj leads to 

oo , 1 \ n 

Ves(E) = V H irj + vHjX—— XHjX—— XH lV . (2.7) 



n =o E-H \ E-H , 

In this form it is obviously identical with the result obtained from "old-fashioned" perturbation 
theory HB 

The effective Hamiltonian for nucleons can also be derived from eq. ( |2,3j ) in another way, 
which we shall refer to as the method of unitary transformation. Introducing new states \x) and 
| if), which are related to \cp) and by the unitary transformation 



Ix) 



one can rewrite eq. ( |2.3[ ) in an equivalent form 

We can decouple the two subspaces \x) and \ip) by choosing U such that the operator WHU is 
diagonal. We adopt the ansatz of Okubo |15| for U. It has the form 



/ (l + ^U)-i/2 -it(i + iit)-iM 

^"^(1+^-1/2 (i +j4 ^t)-i/2 ; ^- iu i 

with the operator A satisfying 

A = XAr]. (2.11) 
The requirement for U^HU to be diagonal leads to the following nonlinear equation for A: 

X (H - [A, H] - AHA) 7] = . (2.12) 

In the case when the interaction Hamiltonian Hj can be treated as a small perturbation, it is 
possible to solve Eq. ( p. 12 ) perturbatively to any order. For instance for the Hamiltonian H 
represented by 

oo 

H = H + Y / H n (2.13) 

n=l 

with the index n denoting the power of the coupling constant, one assumes the operator A to 
be of the form 

oo 

A = Y,A n ■ (2.14) 

71=1 

The solution of eq. ( |2.12| ) to order n is then given by 

^ f n—l n— 1 n— 2 71— J— 1 1 

A n = —A <H n + J2 H i A n-i ~ Y. K-iHi - V £ AiHjA^j \ 7] . (2.15) 

£ ~ Ho 1 7=1 7=1 7 = 1 j = l J 
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Here, we denote the free-particle energy of the state \rf) by £. One can see from eq. ( 2.15 ) that 
it is possible to find A n for every n recursively, starting from A\. The effective Hamiltonian, 
which operates solely in the subspace \x), is given by 

iJ eff = 77(1 + A^A)- 1 / 2 (h + A*H + HA + A* HA) (1 + A^A)- 1 / 2 ^ , (2.16) 



as it follows from eq. (|2^). Expanding (1 + A^A) 1 / 2 and using eqs. ( |2.13| ), Q2.14D , ( |2.15| ) one 



can obtain the effective Hamiltonian to any order in the coupling constant. 

Several modifications are necessary by applying the formalism described above to effective 
Lagrangians (Hamiltonians) and in particular to chiral invariant Lagrangians. First, the expan- 
sion in powers of a coupling constant must be replaced by the expansion in powers of small 
momenta. For doing that, power counting rules are necessary. Furthermore, one expects the 
operator XAr] to consist of an infinite number of terms to any order of Q caused by infinite 
number of vertices in the Hamiltonian. We now show how these problems can be solved. 

3 Application to chiral invariant Hamiltonians 

We want first to recall the structure of the most general chiral invariant Hamilton density for 
pions and nucleons, 

H = Hq+Hi . (3.1) 

We shall treat the nucleons nonrelativistically as it also has been done in [H f|] [f§] • The purely 
nucleonic part of Hq is nothing but the kinetic energy 

n N0 = -N^N , (3.2) 
2m 

where m denotes the nucleon mass. The mass term disappears from the Hamiltonian when one 
separates the basic time-dependence of nucleon fields (in the rest-frame) via 

N = e- tmt N . (3.3) 

In higher orders in small momenta, which we shall not treat here, relativistic corrections to 
eq. ( [3.2D must be taken into account and the free nucleonic Hamilton density looks like 

H N0 = H N { ? + n N f + Hn^ + ... (3.4) 

with Hjvo gi yen by eq. (j3~2j). The upper indices denote the number of derivatives acting on 
N. These terms stem from the expansion of the relativistic kinetic energy in powers of l/2m. 
A path integral approach to derive such relativistic corrections explicitly from the underlying 
relativistic theory in the one-baryon sector is spelled out in |l8}| . The free Hamilton density for 
pion fields is given by 

n^ = \* 2 + \{Vnf + \ m W , (3.5) 

where the denotes the time derivative and m n the pion mass. We split the interaction 
Hamilton density into three parts: 

Hi = H n7T + Hnn + ■ (3.6) 

The first piece describes self-interactions of pions and contains an even number of derivatives 
and pion field operators and any number of m^-factors. The terms with two derivatives (or one 
m 2 ) have at least four pion fields. 
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The second piece in eq. (|3.6| ) consists of four or more nucleon fields and any number of 
derivatives. The terms of 7i n N have any number of pion fields and either two nucleon fields and 
at least one derivative (or m 2 ) or four or more nucleon fields and any number of derivatives and 
factors of . The structure of the effective Hamiltonian is now clarified and the explicit form 
we use is given in Appendix C. Because of baryon number conservation and the absence of anti- 
nucleons, the subspaces of the Fock space with different number of nucleons are automatically 
decoupled. That is why we define the state \4>) as consisting of N nucleons only. The state \tp) 
contains then N nucleons and at least one pion. To avoid difficulties with an infinite number of 
vertices in the Hamiltonian, we shall split the subspace into a series of orthonormal subspaces 
with definite numbers of pions 

\^) = \^) + \^) + \^) + ... . (3.7) 
The indices denote the number of pions. Analogously we introduce the corresponding projectors 

A = A 1 + A 2 + A 3 + . . . (3.8) 

with the properties 

X i X j = 6 ij X i (3.9) 
Using eq. ( |3.9| ) we can rewrite eq. ( |2.12| ) in the form of an infinite system of equations 

A 1 (H - [A, H] - AHA) rj = 
A 2 (H - [A, H] - AHA) rj = 

(3.10) 

X 1 (H - [A, H] - AHA) rj = 



Because of the property eq. (2.11) of the operator A each single equation in eqs. ( |3.10| ) can be 
expressed as 

oo oo 

A^ff/ry + XH^Ar] + ^H^Ar] - XArjHw - KA^H^ - ^ yAqH^Arj = (3.11) 



To solve the system of eqs. ( 3.10 ), we make use of the usual philosophy of effective theories. We 
are interested only in low-energy processes. When the renormalization scale for loops is chosen 
to be of the same order as the low-energy scale Q, it is the only quantity with dimensionality 
of energy apart from the coupling constants. The low-energy matrix elements can then be 
classified by powers of Q by use of simple dimensional analysis. This immediately raises the 



following question: The projection of the Hamiltonian in eq. ( [2.16 ) involves (in principle) an 



integration over all momenta, whereas the matrix elements of the operator A, which can be found 
from eq. ( |3,10| ), are obtained from a low-energy expansion. Does this amount to a consistent 
procedure? The answer to this is yes, because of the renormalization, after which all typical 
momenta in virtual processes are bounded and of the order of Q. 

Now we make an ansatz for the operator X l Arj, which will be justified at the end of this section: 
we assume that it consists of equal number of vertices from Hj and energy denominators. Each 
of these denominators contributes a factor of 1/Q- It is now easy to calculate the power ua of 
Q of the matrix element {tjj l \A\(j)) using the usual topological identities, 

(il> i \A\4>)=Q v *r(Q/A x ,n/A x ,g i ) , (3.12) 
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with \i the regularization scale (since regularization is needed to render certain loop integrals 
finite) and gi collectively denoting all coupling constants (low-energy constants). Each integra- 
tion over virtual momenta brings three powers and each derivative one power of Q, respectively. 
The pionic phase-space factors l/^2u>, the energy denominators and ^-functions give each a 
factor of Q~ 1//2 , Q^ 1 and Q~ 3 , in order. For L loops, C separately connected pieces, I p internal 
pion lines and Vi vertices of type i with di derivatives one finds 

v A = 3L + ]T Vidi -I p -^2Vi-3(C- 1) (3.13) 

i i 

According to our assumption J2i Vi stands for the number of energy denominators. Here and in 
what follows we do not count the overall (^-function and the phase-space factors of the external 
pions. Now we use the general identities 

L = C + I p + I n -J2Vi , (3.14) 

i 

2I n + E n = Y J V i n i , (3.15) 

i 

and 

2I p + E p = J2ViPi , (3.16) 

i 

where E n = 2N is the number of external nucleon lines, E p the number of external pion lines, I n 
the number of internal nucleon lines and rii (pi) the number of nucleon (pion) fields at a vertex 
of type i, to cast eq. ( |3.13| ) into the form 

v A = 3 - 3N - E p + Vi «i , (3-17) 

i 

with 

3 

«i = di + -m +pi-4 . (3.18) 

Again we point out that this result takes its form due to the ansatz that the projected operators 
\ l Ar] with a given number of nucleons and pions consist of an equal number of vertices and 
energy denominators. Also, it should be mentioned that we use the number of external pion 
lines to express the counting index va instead of the number of loops and of separately connected 
pieces as it has been done in [|lj]-[|5[]. This is more natural in the projection formalism employed 
here. Let us take a closer look at vertex dimension Ki, which is related to the canonical field 
dimension. It is well known that all interactions can be classified with respect to Ki. Those 
with m < are called relevant, with Ki = marginal and with Ki > irrelevant. The last ones 
are "harmless" and can be well treated within low-energy effective field theories. In contrast 
to these, the relevant and marginal interactions lead to complications with the power counting. 
This becomes immediately clear when one takes a look at eq. ( p. 17 ): the number of possible 



diagrams for a given process and to a given order is no more bounded. Furthermore, for relevant 
interactions the number 1/4 is even not bounded from below. That is why the perturbative 
treatment is not possible in this case. For more details about the role of such interactions 



in effective field theories see [19|. Chiral symmetry does not allow any relevant or marginal 
interactions in the interaction Hamiltonian Hj. The minimum possible value of Ki for vertices 
in Hi is one. The corresponding vertex contains two nucleons, one pion and one derivative. 
In general, for vertices with two nucleons and pi pions one can easily find that Ki > p%. The 
interactions with pions only have the vertex dimension K\ > 2 with k\ = 2 corresponding to 
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the interaction with four pion fields and two derivatives. Again, for those vertices with an even 
number of pions and without nucleons it follows, that ftj > max(2, — 2 + pj). Using eq. ( 3.17| ) 
one can determine the minimum value of 1/4 for \ 4k+t Arj with k positive integer or zero and 
i = 0,1,2,3: 

min(^ A ) = 3 - 3N - 2k . (3.19) 

All diagrams corresponding to such operators consist of maximum possible number of pure pionic 
vertices with pi pion fields and Ki = — 2 + pi and any number of interactions with two nucleons, 
Pi pions and Ki = p% without intermediate pion states. In the case i = 3 one intermediate pion 
state may also occur. 

We define the order of the operator X ik+t Ar] by the number I (I = 0, 1, 2, 3, . . .) which satisfies 



v A = 3 - 3N - 2k + I 



(3.20) 



and introduce the following notation: \ ik+l Air]. Let us now find the power v of Q for every 
term in eq. ( 3.11| ). By H K we denote a vertex from the interaction Hamiltonian Hi, with the 
index k given by eq. (|3.1S| ). Details are relegated to appendix A. It can be seen from eqs. ( |A.7| ), 
flAlCl ), flO|), QA.15D and flP6| ) that the minimum possible value of v for eq. ( |3.11| ) is given 
by 

min(z/) = 4 - 3N - 2k . (3.21) 

We now show how to solve the system of eqs. ( [3.10 ) perturbatively. For that, we define the 
order of eq. ( |3.11 ) to be r > given by 



v - (4 - 3N - 2k) 



(3.22) 



In appendix B we check what kind of operators X 4k+l Air] enter the eq. ( 3.1 1] ) at each order r. 
That equation can be expressed as 

E(\ Ak+i )\ Ak+i A l=r rj = -X 4k+i {H - [A, H] - AHA} rj . (3.23) 

The free energy is always denoted by E. In this case it is the pionic part only. In the curly 
brackets we have written in symbolic form all remaining terms of eq. ( 3.11 ). For i = 1,2 they 
contain only operators of the type \ Ak+t Aii] with 



Ak + i < 4k + i for I 



or I < r 



(3.24) 



As an example we reg ard \ 4k+i HX ik+i Air]. For 4k + i < Ak + i it follows from eqs. (Q, flB~3p 
that / < r. For 4k + i = 4k + i, eq. (|R4|) leads to I < r — 2. Finally for 4k + i > 4k + i 
eqs. (B.5)-(B.8) require I < r — 2. Eq. ( |3.24| ) is also true for i = 3 apart from the single term 
_ A 4fc+3# iA 4(fc+i)^ =r?? _ Thig ig due tQ eq _ In the 

case i = 0, only the operators \ 4k+i Air], 

restricted by the conditions 



0, k < k for / 



or I < r 



(3.25) 



can enter the right hand side of eq. ( |3.23j ). In all cases the value of k is bounded from above by 
the inequality 

k < k + -(r -1 + 8) . (3.26) 
6 



as can be defered from the second inequality in flB,8| ) and the equality in ( B.16| ). Now it is 
clear how to deal with the system of equations ( |3.10| ). At each order r one starts by solving 
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the equations with i = and with increasing number k to define all operators X 4k Ai =r rj. This 
requires knowledge of operators X^Airj of the same order I = r with k < k and a finite number of 
operators in lower orders I < r. The next step is to solve the remaining equations with i = 1, 2, 3 
with increasing number Ak + i starting from Ak + i = 1. After that one can go to the next order 
r + 1. The number of equations to be solved in each order can be estimated by use of eq. (ET26D 
and the inequalities of appendix B. 

To justify our ansatz about the structure of the operator A, consider the starting equations 
at order r = 0, given by 

E{X 4 )X 4 Aov = X 4 H 2 t] , (3.27) 
E(X 1 )X 1 A r] = X l Hi-q , (3.28) 

where E(X l ) corresponds to the kinetic energy of the pions in the state A\ From eq. ( |3~2^ ) 
one can see that the unknown operator X ik+l Ai =r r] has the structure which we assumed at the 
beginning of this section, if and only if the already known operators XAr/ entering the right 
hand side of this equation have precisely this form. Therefore, to proof our ansatz recursively 
for all operators XA-q it is sufficient to see that it holds for the corresponding starting operators 
in eqs. (|3.27 ) and ( 3.28 ), which is obviously the case. 



4 Two— and three— nucleon forces using the method of unitary 
transformation 

We shall now practically apply the formalism described in the last section and derive an effective 
Hamiltonian for nucleons in leading and next-to-leading orders. As a starting point we use the 
effective chiral invariant Hamiltonian for nucleons and pions, which was specified at the beginning 
of the last section. Its explicit form is given in appendix |C| 

4.1 Two and many nucleon potential in the projection formalism 

First of all we look at the eq. ( |2.16| ) and define the order of all terms on its right hand side. 
Introducing the operators 

L t = ^Ajx^+^A^r, , (4.1) 
and its chiral power (the power of Q) 

v t = 3 - 37V + v t (4.2) 

with 

v t = Ak t + 2i t + l t + l' t , (4.3) 
evaluated via eqs. ( |3.17 ) and ( |A.1| ) one gets contributions of the following types. In the following 



v denotes the corresponding power of Q. 



' [n^ 



t 

V 



3 - 3N + J2 h (4.4) 

t 

2. [n^]^[n^ 



9 



1/ = A -3N + K + J2"t + J2 i> s ( 4 - 5 ) 

t s 
t s 

and [l[L t ]rjH K X Ak ^A lmV [l[L s ] 

t s 

v = 4 - 3N + k + 2k m + i m + l m + vt + ^ Ds ( 46 ) 

t s 

4. [ I] Lt] nAlx^+^H^+^A^ [ J] L s ] 

t s 

v = 4 — 3N + K + l m + l n + 2k m + 2k n + i m + i n + j> t + u a . (4.7) 

( s 

In case of the unit operator in eq. ( |2,16| ) one has to drop the 2?'s. The effective potential can be 
easily read off from the effective Hamiltonian via 

V eS = H eS -H . (4.8) 

From now, we denote the contribution from the free Hamiltonian Hq by the pertinent nucleonic 
and pionic free energies, £ and u, respectively. Using now eqs. ( ^.16|) and (fO|) - (|4.7|) one can 
write the leading term of order v = 6 — 3iV of the A^-body potential as: 

V^- 3N) =v(H 2 + AX 1 ^ + H^Aq + A^coAo) r, . (4.9) 

Note that in the last term only the pionic free energy u> contributes in lowest order. As noted in 
|H and |§, there are is contribution to the potential to order 7 — 3 A?" because of parity invariance. 
Another way to see this is the fact that there is no nonvanishing operator X l H 2 rj. At next order, 
8 — 3N, one gets a more complicated expression for the potential: 

VS~ 3N) = r ] (H 4 + Al\ 1 H 1 + Alx 1 H 3 + H 3 X 1 A + H 1 X 1 A 2 



+ 4a 1 F 2 A 1 A - -AJaMqt^ - -HiqA^M 

+ Alx l £A Q - ^A^AqS - -SA^Aq + A^uAo + A ] X 1 uoA 2 

+ a\x 1 H 1 X 2 Aq + J^X 2 HiX 1 A q (4.10) 

- ^A^AovAlX 1 ^ - ^Alx'Aov^X'Ao - ^X 1 H^X 1 A 

- ^H^A^A^Ao - ^Alx l A or] Alx l ojA - ^A^uA^A^Aq 

+ A ] Q X 2 H 2 + H 2 X 2 A + 4 A Vi + w 2 )4)) V ■ 

Note that the S's denote the nucleonic free energies related to the accompanying projection 
operators (A or 77). The operators X 1 Aqtj, X 2 Aq7], X 1 A 2 t] can be evaluated along the lines 
described in the last section. This leads to: 

A 1 Aw = -—H lV , (4.11) 

to 

\2 \2 

X 2 A oV = ■ H 2V ■ H^Aqtj, (4.12) 

U)\ + UJ 2 U)\ + U) 2 
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A 1 A 1 A 1 A 1 
X 1 A 2 r J = #37? H^Aorj H 2 X 1 A 7 ] + —A r ] H 1 \ 1 A r ] (4.13) 

UJ UJ UJ UJ 

A 1 A 1 A 1 

+ —A oV H 2 7] - £—A 7] + — A oV £ . 

UJ UJ UJ 

Putting eqs. (4.11) - ( [4.13| ) into eqs. ( f4.9[ ) and (4.10) and performing straightforward algebraic 
manipulations, one obtains the potential expressed as 

Vt m) = V^-H^HX (4.14) 

\ CO ' 

<~ m = r,(tf 4 -tfi-tf3-#3-tfl 

v UJ UJ 

A 1 A 1 A 1 A 2 A 2 A 1 
+ H\ — H 2 — Hi + Hi — Hi H 2 + H 2 Hi — Hi 

UJ UJ UJ UJl + UJ 2 UJl + UJ 2 UJ 

- H 2 A ' H 2 (4.15) 

UJl + UJ2 

— Hi — Hi ; Hi — Hi + -Hi—-tHi7]Hi — Hi + -Hi — Hir]Hi——^Hi 

uj uji + uj 2 uj 2 [ujY uj I uj [ujy 



1 A 1 1 A 1 

— ~H\ - — —^H^TlHo — — HoTlHi - — rrrH) 

2 (uj) 2 1 2 ' (uj) 2 

[uj) z 2 (uj) z 2 (uj) z > 

Consider first the leading order potential, eq. ( 4.14j ). It consists of two terms, the one-pion 
exchange ~ Hi(X 1 /uj)Hi and the four-nucleon contact interactions subsumed in H 2 . This 
potential obviously agrees with the one obtained in time-dependent perturbation theory. More 
interesting is the first correction given in eq. ( gig ). Thedia grams corresponding to the various 
terms are shown in figs. 1-4 and 6. The first line refers to the graphs of fig. 1, the second to 1, 
2, 3 in fig. 2 and 1, 2 in fig. 6 and the third to graph 4 of fig. 2. We should mention that all 
graphs containing vertex corrections with exactly one 7T7riViV-vertex, which are also contained 
in the second line, give no contributions, because only odd functions of the loop momentum 
enter the corresponding integrals. The first term in the fourth line subsumes graphs 5 to 8 of 
fig. 2 plus the irreducible self-energy diagrams depicted in fig. 3. The next two terms in the 
fourth line refer to graphs 9 and 10 plus the "reducible" self-energy diagrams of fig. 4. Such 
"reducible" diagrams are typical for the method of unitary transformation and do not occur 
in old-fashioned time-ordered perturbation theory. They should not be confused with truly 
reducible diagrams, one example being shown in fig. 5a. In that figure, the horizontal dashed 
lines represent the states whose free energy enters the pertinent energy denominators. In old- 
fashioned time-ordered perturbation theory, such reducible diagrams are generated by iterating 
the potential in a Lippmann-Schwinger equation, with the potential being defined to consist only 
of truly irreducible diagrams. In contrast to the really reducible graphs like the one in fig. 5a, 
the ones resulting by applying the projection formalism do not contain the energy denominators 
corresponding to the propagation of nucleons only. In the same notation as used for fig. 5a, one 
can e.g. express diagram 9 of fig. 2 as a sum of two graphs as depicted in fig. 5b, c. All diagrams 
involving four-nucleon contact interactions, shown in fig. 6, follow from the fifth line and, as 
already noted above, from the first term in the second line of eq. (4.15). We remark that in the 
c. m. system the three terms in the last line add up to zero. We have nevertheless made them 
explicit here since in time-ordered perturbation theory, these terms are treated differently and 
lead to the recoil correction, i.e. the explicit energy-dependence, of the two-nucleon potential. 
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4.2 Two— nucleon potential: Expressions and discussion 

As already noted above, the leading part of two nucleon potential is given by just one pion 



exchange with both vertices coming from eq. ( C.l ) and the nucleons treated as static particles 



plus contact interactions corresponding to eq. (C.3). This result holds in both time-dependent 



perturbation theory and projection formalism. For easier comparison, in this paper we use the 
same notation and definitions as it has been done in |Q] . In particular, the spin and isospin 
matrices a and t satisfy the following relations: 

Utj = jSij + ^ie ijk tk , (4-16) 
<Tj(Tj = 5ij + ieijkcrk ■ (4-17) 

The initial (final) momentum of nucleon in c. m. system is denoted by p (p')- The transferred 
and average momenta are given by q = p — p' and k = (p + p')/2, respectively. The leading 
order potential is then given by 

yV> = -( 2 -p-) 2 tl . t 2 fr;pV?> + C S + C T B X ■ B 2 (4.18) 
V F n J q 1 + ml 

The first corrections to this result appear at order two. We now enumerate all such corrections. 
We start with the ones obtained using old-fashioned time-dependent perturbation theory HQ. 
First, one seemingly obtains the contributions of tree diagrams 1 and 2 in fig. 1 with one vertex 
coming from eq. ( |C.5| ) and still using the static approximation for the nucleons, 

<U = • t 2 0^tp> (A ig -> + A 2 f) . (4.19) 

The constants A\ and A 2 are given by 

Art \ , . { .1 Ar, 



A! 



A[ + ^) , A 2 = -4 [a\ - ^) . (4.20) 



We note that these expressions differ from the ones given in appendix B of H. Since k and q 
are defined with a different prefactor as sum and difference of p and p', there can not be the 
same prefactor relating the (^1,^2) to the (A^A^). A few more remarks are in order: If one 
decomposes q as q + ml — m|, one sees that the ^-contribution can be completely absorbed 
in coupling constant and contact term redefinition. However, as discussed in appendix C, the 
term ( |C.5| ) appears in the Hamiltonian with a fixed coefficient ~ 1/m representing relativistic 
corrections and is not relevant to the order we are working if the nucleon mass is counted via 
eq. (IA.14D . For this reason the contribution ( [4.19| ) to the two-nucleon potential can and should 



be completely omitted. 

The contributions from tree graphs with various contact interactions from eq. ( C.4 ) can be 
expressed as 

Vyltree = C x q 2 + C 2 f + {C z q 2 + C^\o 1 ■ ff 2 ) + iC$^^ ■ {q * k) 

+ C 6 (g • <?i)(g • a 2 ) + C 7 {k ■ a x ){% ■ a 2 ) (4.21) 
where the Cj's are given by 

Ci = c[-c> + ^, 
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C 2 = AC' l -AC' z -2C' 2 , 
r" 

3 — Og + — U 14 , 

C 4 = 4C 9 - 2C( 2 + 4C( 4 , (4.22) 
C 5 = 2C' 5 -4C' 4 -2C' 6l 



£10 + Wi _ c / 5 
2 2 u 

C 7 = 4C^ + 4C' 8 - 2C[ - 2C' n + iC[ 3 . 



C& - Cj + C' 8 + + — - , 



Note again that the contributions ~ C5 , ~ C( 1 and ~ C 14 and also the signs in many places are 
different from what is given in Q . Further corrections arise from irreducible one-loop diagrams 
1-8 in fig. 2: 



/r2 

^(2) _ *9A 



V 2tt. 1-loon. irr. ~~ mA / tajl 1 !'^ 



27 r,l-loop,irr. p 4 J ( 2 7r) 3 ' 



1 /■ d 3 Z (w+-w_) 2 1 

' -tl-t 2 — (4.23) 



2F^ 7 (2vr) 3 cj + + w_ 

^ /■ d 3 / 1 j ( 3 | 8ti • t 2 \ y . 2 \ 2 



+ 4 ( j + • [?x r|)(ai • [gx f]) 



where 

u± = yj {q± I) 2 + Ami • ( 4 - 24 ) 

We have also calculated the contributions from various irreducible one-loop diagrams 1-8 in 
fig. 3 and 1-2 in fig. 6, which involve self-energy insertions and vertex corrections. These have 
not been considered in HQ. The corresponding contributions are given by 



+ ±)f(a l .q)(a 2 .q)\ , (4.25) 



where 



and by 



q 2 + ml , wj = V f + ml , (4.26) 



41l-loo P ,ir, = "ffs J ^3^{3r 2 + 4(t 1 .t 2 )(r.a 1 )(r.a 2 )} (4.27) 

I tPy^ \Ufti-t 2 + 2{l-a l ){T-a 2 )(?,-2t 1 -t 2 ) 
F 2 J (2vr) 3 uf { 

-f{ai •o f 2 )(3-4ti -t 2 )} • 

(2) -.2 
We remark an unpleasant feature of l-loop irr ■ With t = —q one can rewrite the term 

proportional to as (ml — i)~ 3 / 2 . This function has a cut starting at t = ml- Physically, 
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this does not make sense since all cuts should be produced by multi-particle intermediate states. 
We come back to this later on. 

Finally, the corrections arising by explicitly keeping the nucleon kinetic energy in one-pion 
exchange tree diagrams with both vertices from leading order Lagrangian represents the energy- 
dependent part of 2N-potential given by 



VP = ^ • t 2 {a x ■ q)(a 2 ■ q)—^ , (4.28) 

(q 2 + mlj 



where E is an initial energy of the two nucleons (full energy). It should be stressed that on- 
shell, this contribution vanishes. This on-shell center-of-mass kinematics is often used as an 
approximation for calculations of NN scattering or few-nucleon forces. However, if one iterates 
this potential in e.g. a Lippmann-Schwinger equation, the energy E can no longer be simply 
related to the momenta k and q . A similar comment applies to using this recoil correction in 
the calculation of few-nucleon forces (for N > 3). 

Let us now look at the corrections arising in the framework of the projection formalism. The 
corrections from tree diagrams with one vertex taken from next-to-leading order Lagrangian 
do not change and are given again by eq. ( 4.2l| ). Remember that we do not consider the 



Hamiltonian H.3 given by eq. (|C.5[) for reasons given above and in appendix C. Consequently, 
there is no contribution ( [4,19| ) to the potential. Apart from the corrections from irreducible 
one-loop graphs 1-8 in fig. 2 given by eq. ( 4.23| ) one obtains a contributions from reducible 
diagrams 9 and 10, which can be expressed by 

V 1-loop, red. " / J^3 [ J ~ 2 *1 " *2 J ^3 i 4 ' 29 ) 

a 2 - [qx f])(ai • [q x f]) \ . 




Summing up the eqs. (|42p and ([4291) we obtain the two-pion exchange contributions to the 



potential within the projection formalism: 

v m _ f dH {? - f) 

^,1-loop - p4 J (2vr)3 tl - t2 W+a ;_( a ; ++W _) 

d 3 l 



t x .t 2 ^ + - U - )2 1 (4.30) 



2F£ J (2vr) 3 uj+ + u- lu+lj 

9a f d31 ^1_+uj+uj- +ul J 0+ + f f 2 _, 2 ^ 
+ 



F4 J (2tt) 3 + u-) \ 2 V 

+ 3(a 2 -[qxl])(a 1 -[qxf]) \ . 



As was first noted in ref. |^0| using a different formalism, the isoscalar spin independent central 
and the isovector spin dependent parts of the two-nucleon potential corresponding to the two- 
pion exchange adds up to zero. It is comforting that we find the same result. Note, that this 
is different from the energy-dependent potential derived using the time-ordered perturbation 
theory. 
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The corrections arising from one-loop reducible diagrams 1-4 in fig. 4 and 3-6 in fig. 6, which 
involve the nucleon self-energy and loop corrections to the four-fermion interactions, are given 
by 

(2) '-.7.1 / " ' , , ^/ 



l7r, 1— loop, red. 



and by 



(2) 



NN, l-loop,red. 



1M 

F 4 



9 a r 

IT 



d 3 l UJl+UJ q f2 

<i *2 — -3-3 1 \ a i ■ Q) \ a 2 ■ q , 



(27T) 



(4.31) 



F 2 



Ct 



d 3 l 1 r -q - i 

___| 3 ; + 4(t 1 -t 2 )(l-a 1 )(l-a 2 )} (4.32) 

(i 3 / 1 f -2 
-—3 — 4Z t x -t 2 + 4(ti • t 2 )(Z • <?i)(/ • cr 2 ) 



-Z (^i-^ 2 )(3-4t 1 -t 2 ) 



Note that 



(2) 



lw, 1— loop, red. 



again contains a term 



1K 3 



Combining eqs. (|4.25) and fl4.31| ), this 



unphysical contribution vanishes. Summing up the corrections eqs. ( 4.25| ), (|4.27 ), correspond- 
ing to irreducible graphs, which are again the same as in the old-fashioned time-dependent 



perturbation theory and those (|4.31| ), ( |4.32[) , corresponding to reducible diagrams, one gets the 
complete result, representing the one-loop contributions, which involve self-energy insertions 
and vertex corrections within projection formalism: 



v, 



(2) 



V, 



Itt, 1— loop 



(2) 

NN, 1-loop 



45i 
F 4 



7^-T3 —3-^" 1 (* ' 9 ) ( ■ * ) (°2 ■ 9 ) + m ■ Q ) (0-2 • I 



+22 (#i-q)(ff 2 -q) 



F 2 



Ct 



d 3 l 1 



(2tt) 2 



ti-t 2 ) (I ■ ■ a 2 ) 



<*2. 



(4.33) 



(4.34) 



The complete potential of second order in the projection formalism is therefore given by the 
expressions flOID , flOcD , ( |433|) and g3|). 

We remark that one can perform the /-integrations in the last two equations and one finds that 



the contribution from eq. (4.33) can be entirely absorbed in coupling constant and wave function 



renormalization whereas eq. (4.34) merely amounts to the renormalization of the lowest order 
four-fermion contact terms. Finally, the crucial difference between the two formalisms consists 
in the treatment of the energy-dependent term eq. (4.28). As already noted above, it does not 
appear in the method of unitary transformation. We note that many of the results derived here 



have already been found in [ 20 [ , where one- and two-pion exchange graphs were calculated by 
means of Feynman diagrams and using dimensional regularization. That potential can and has 
been applied for scattering processes but can not be used to solve the bound-state problem. 



4.3 The three— nucleon potential 

Let us consider the leading part of the three-nucleon (3N)— potential. For this we take a closer 
look at eq. ( |4.15 ). The diagram 1 in fig. 7 represents all terms in the second line of this equation, 
which involve the nonlinear ttttNN coupling. These contributions are identical in both time- 
dependent perturbation theory and in projection formalism and, as it was noted in ref. Q, sum 
up to zero, when all time-orderings are taken into account. For three-nucleon forces mediated 
by two-pion exchange with only linear ttNN couplings one gets different results by using the 
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two approaches. The irreducible graphs 1-4 in fig. 8 represent the result of old-fashioned time- 
dependent perturbation theory. Apart from them, the reducible graphs 5 and 6 in fig. 8 occur 
using our approach. An exact cancellation between the contributions from those reducible and 
irreducible diagrams was recently pointed out in [[0]] for the case of an expansion in the pion- 
nucleon coupling constant and adopting the nonrelativistic approximation for nucleons. At 
leading order, these two very different expansions lead to the same set of diagrams. We remark, 
however, that only for the method employed here one has a consistent scheme of treating the 
next-to-leading order corrections. Note furthermore that fig. 8 differs from the analogous one 
given in [T^j by the topology of diagram 6. The graphs e and g shown in fig. 3 of |ll| lead in 
general to different prefactors multiplying the energy denominators. To understand why these 
cancellations appear, let us take a closer look at the terms in the fourth line of eq. Q4.15 ), We 



just rewrite them schematically by pulling out the common factor M, representing the spin, 
isospin and momentum structure, which is obviously the same for all of these terms. For the 
irreducible diagrams we then have 



+ : r + 



Wl(wi + UJ 2 )UJ2 + OJ2) (UJ\ + 



M = -^±^M . (4.35) 



LOtiO. 



1^2 



The contribution from the reducible diagrams can be expressed as 



1 



Uj\u2 OJ1U2 



M = ^±^M . (4.36) 



The cancellation is now evident. The same sort of cancellation can be observed for diagrams 2, 
3, 4 in fig. 7 and for 1-loop graphs 1-6 in fig. 6 with the contact vertex (l/2)Cs(N^ N)(N^ N) 
contributing to the two-nucleon potential, which correspond to the first term in the second line 
and the terms in the fifth line of eq. ( |4.15| ). We conclude that there is no three-nucleon force 
at the order u = 1. This agrees with the finding of Weinberg using time-ordered perturbation 
theory JjJ]. In that case, the mechanism is different. Although the leading three-nucleon force 
corresponding to diagrams 2 in fig. 7 and 1-4 in fig. 8 does not vanish in that case, it was noted 
in |l[ that these contributions cancel exactly with the energy-dependent part of the two-nucleon 
potential when the latter is iterated in the Lippmann-Schwinger equation. That is the reason 
why one can describe systems of three or more nucleons using the energy-independent part of 
the 2N-potential and without explicit three- and many-nucleon forces to the order in the chiral 
expansion also considered here. However, such type of cancellations does not help to remove 
the problems in the two-nucleon sector due to the explicit energy-dependence as noted before 
(i.e. that the wavefunctions are only orthonormal to the order one is working). This is the main 
difference to the projection approach we are using. 



5 Summary and conclusions 

In this paper, we have presented a novel approach to the problem of deriving the forces between 
few (two, three, . . .) nucleons from effective chiral Lagrangians. For that, we use the well-known 
projection formalism dating back to Okubo, Fukuda, Sawada and Taketani [15|[16|. It allows 
to uniquely separate the effects from the operators and the wave functions. For the case at 
hand, we first had to modify the power counting rules since in the projection formalism one 
decomposes the Hilbert space of nucleons and pions into subspaces with definite nucleon and 
pion numbers. While in old-fashioned time-ordered perturbation theory the resulting wave 
functions are only orthonormal to a certain order in the chiral expansion, this problem does not 
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occur in the projection formalism. Furthermore, in the previous calculations based on time- 
ordered perturbation theory, the two-nucleon potentials turn out to be energy-dependent, which 
is a severe complication for applying these in systems with three or more nucleons (although to 
leading order, these recoil corrections are cancelled by certain N-body interactions). Here, we 
have given the explicit expressions of the two-nucleon potential to two orders beyond the lowest 
order tree graphs (comprising the leading one-pion exchange and four-nucleon contact terms). 
In particular, we have also calculated self-energy type corrections and one-loop corrections 
to these four-fermion operators not considered before. The most salient features of the so 
constructed two-nucleon potential is its energy-independence and the orthonormality of the 
corresponding wave functions. Furthermore, we have considered the leading contributions to 
the three-nucleon potential. As in time-ordered perturbation theory, we find that these sum up 
to zero. However, the mechanism for the cancellation between some of the graphs is distinctively 
different in the projection formalism since it is not related to an intricate cancellation of terms 
generated by the iteration of the energy-dependent two-nucleon potential in old-fashioned time- 
ordered perturbation theory. Rather, in the approach here, these cancellations can be traced 
back to the appearance of "reducible" graphs whose precise meaning is explained in section 4.1. 
These diagrams are in fact responsible for the orthonormality of the wave functions and are thus 
sometimes called "wave function re-orthonormalization" graphs. 

Clearly, the results presented here constitute only the first step in a bigger program. We are 
presently working on a regularization scheme, which allows us to consistently treat the potential 
and the Lippmann-Schwinger equation to generate the bound state (for a recent discussion with 
many references to other approaches, see e.g. Beane et al. |^l|). In the pioneering work by van 
Kolck et al. ||], this was not done but rather a Gaussian cut-off for the momenta q and k was 
introduced irrespective of the power counting. Furthermore, the A (1232) degrees of freedom 
have to be incorporated in a manner analogous to what was done for the single nucleon case 
by Hemmert et al. |22| and finally, numerical results have to be obtained for the two- and 
three-nucleon systems. We believe that within the framework outlined here, one can find out to 
what extent chiral symmetry plays a role for the nuclear forces (the role of one- and two-pion 
exchanges for the higher NN partial waves was recently elucidated in [^3] |2(]]). Furthermore, this 
approach can also be used to calculate meson-exchange currents in a manner which is completely 
consistent with the so-constructed nucleon potential. 
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A Chiral counting for the projected equations 



In this appendix, we work out the chiral power of the different terms appearing in eq. ( 3.1l| ). 



Any one of these can be characterized by a counting index i.e. is proportional to Q v . This 
index should not be confused with the one giving the chiral dimension of the matrix-elements 
discussed in the text. The free part of the Hamiltonian Hq gives the energy, which we 
denote throughout by E. Remember that H K denotes a term from the interaction Hamiltionian 
Hi with the vertex dimension n. 

One can determine the chiral power v for all terms entering eq. ( |3. 11 ) by using of the cor- 



responding result ( 3.17] ) for va, which, however, should be modified to take into account the 
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correct number of energy denominators in each case, because of the presence of the operators Hq 
and Hj. Furthermore, one uses eqs. ( 3.17] ), ( |3.20| ) to express the sum J2i Vi K t over an vertices 
entering the operator X Ak+l Air) in terms of k, i, and I as follows: 



=2k + i + l . 

i 

For example, for X 4k+l H K \ 4q+ ^ A^ one has 

v = 3 - 3N + 1 - (4k + i) + k + y i K i =4-3N-4k-i + 2q + j + l + K 



(A.l) 



(A.2) 



where (4k + i) corresponds to the number of external pion lines and we have added +1 because 
there is one energy denominator less than it was counted in eq. ( 3.17| ). Proceeding analogously 
one can derive a corresponding chiral power for all other terms in eq. (ETIlD , which we shall use 
in what follows. 

Now we take a look at the possible vertices in Hi which are explained in section 3 with p pions 
in order to obtain few useful inequalities for determining the corresponding vertex dimension k: 



K > 2, for p = , 
K > p, for p = 1, 2, 3 



(A.3) 
(A.4) 



These follow from ( |3.18j ). For p = n must be larger equal 4 and for p = 1,2,3 n must be 
greater equal 2 and there must be at least one derivative. Further 



n> 1 , 

K > -2 + p . 



(A.5) 
(A.6) 



Whereas ( |A.5| ) is generally valid, (A.6) follows from simple inspection of eq. ( |3.18j ): for n = 
there are at least two derivatives and for n > 2 it is immediately apparent. 
We find the following terms with the corresponding indices v: 

• X Ak+i H K r] 

u = 4-3N-4k-i + K. (A.7) 



This follows from eq. ( 3.17 ) noting that there is one energy denominator less. It follows 
from (Q, ([AID , fljj ), that 



(A- 



v > 4 - 3N, when k = 
v > 2 - 3N, when k> 

\ 4k+i H K \ A i + iAiri 

v = 4-3N-4k-i + 2q + j + n + l . 
The following cases are to be considered: 

1. q < k - 1 

In this case the number of pions p > 4k + i — 4g — j at the vertex. Therefore the 
inequality ( |A.6| ) takes the form 

k > -2 + 4fc + i - 4q - j . (A.9) 

It then follows from eq. (A.8), that 

v>4-3N-2k . (A.10) 
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2. q = k 
Then 

p>\j-i\ . (A.ll) 

We apply eq. ( A.4) to this case and using the (obvious) inequality j — i > — \j — i\ 



one can justify eq. flA.10 ). 
3. q > k + 1 

One can easily prove eq. ( [A. 10 ) for this case using eq. ( |A.5| ) and by noting, that 
j-i> -3. 



E(X 4k + i )\ 4k + i AiT] 



v>A-3N-2k + l>4-3N-2k . 



(A.12) 



This follows immediately from eq. ( 3.20 ) Here the first inequality becomes the exact equal- 
ity when only the pion kinetic energy is taken into account in E(X 4k+1 ). 



\ Ak+i A ir iH K ri 

v = 4-3N-2k + l + K . (A.13) 

There is no nonvanishing operator r]H K r] with k < 0. That is why the inequality ( A.10| ) is 
again valid. 

We shall count the nucleon mass in the same way as it has been done by Weinberg 0: 



m A1 

A. 



(A.14) 



The nucleon kinetic energy must then be counted as Q 3 and we get the following result: 

v>6-3N-2k + l>6-3N-2k . (A.15) 

. X^A^H.X^A^ 

v = 4 - 3iV - 2k + 2q + j + h + l 2 + k > 6 - 3iV - 2k . (A.16) 

B Operators at order r in Eq. (3.11) 



In this appendix, we work out in detail which operators actually appear in eq. ( p.ll ) at order 
r, which is defined in eq. (gjH ). We find 



X Ak+i H R X A i + iAiri 

One obtaines from ( [A.8Q and ( p.22 ) the following identity: 

/ = r + 2{k — q)+i— j — k 
Let us first consider the case i ^ 0: 
1. q < k - 1 



(B.l) 



Using eqs. (|A.9[ ), (B.l) one gets 

I <r -2k + 2q + 2<r . 



(B.2) 
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2. q = k 

(a) j < i 



It follows from eqs. ([A. lift , (|B3|) , that 

Z < r- — | j?" — « | — |— z — j/' < r" . 



(B.3) 



(b) j = i 



In this case eq. (B.l) takes the form 

I = r — k < r — 2 . 



(B.4) 



due to the inequalities (|A.3[ ), (|A.4[ ) and (|A.6[ ). 
(c) j > * + 1 

The first inequality of eq. flB.3|) leads to 



3. q = k + l 

(a) t = 1,2 



/ < r - 2j + 2i < r - 2 . 



(B.5) 



It can be seen from the inequalities ( |A.4[ ), ( A.6 ) that in this case k > 2. Using 
the inequality i — j < 2 we get from eq. ( p.l| ) 



I < r - 2 . 

(b) % = 3 

The general inequality ( |A.5| ) leads immediately to 

i. j = : I < r , 

ii. j > 1 : I < r — 1 



(B.6) 



(B.7) 



4. q>k + 2 

Using eq. (|A.6| ) with the number of pions p given by p = Aq + j — 4/c — i we obtain 
from eq. ( |B.1| ) 

Z < r - 6(<? - k) + 2 + 2i - 2j < r - 6(q - k) + 8 < r - 4 . (B.8) 

The case i = can be considered analogously: 

1. q < k - 2 ^ Z < r- 2 , 

2. q = k - 1 . 
(a) j > 1 



Putting p = 4 — j in eq. ( |A.4| ) we see from eq. (B.l), that 

Z < r - 2 . 

(b) j = 

In this case we have k > 2, as it follows from eq. (|A.6|). We obtain 



(B.9) 



I < r 



(B.10) 
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3. q = k 

It follows from eqs. ([A .3]) and (|A.4|), that 



/ < r - 2 . 

4. q>k + l 

We apply the inequality (|A,£| ) with p = 4q + j — 4k to eq. flB.l|) and obtain 



(B.ll) 



\ Ak+i Air]H K r] 

Here we have used eq. ( 
\ 4k+i Air]E(r]) 



l<r-6(q-k)+2<r-4 . 



Kr . 



I = r — k < r — 2 . 



I < r - 2 . 



h + l2 = r — k — 2q — j<r — 2 



(B.12) 
(B.13) 
(B.14) 

(B.15) 
(B.16) 



C Explicit form of the Hamiltonian 

In this appendix we give the explicit form of interaction Hamiltonian density, which we use in 
this paper. It is based on the effective Lagrangian given in [Q] with an important modification 
to be discussed below. 



Hi 



—Nka ■ -VttN 
2 

NH • (7T X 77)iV 



F 2 

+ ^C T (iV f aiV) • (iVW) + ^C s (iV f iv) (n*N 
C[ [(A^ViV) 2 + (VN^N) 2 } + C^N^VN) ■ (VN^N) 
+ C' 3 (N^N) \N^V 2 N + ^N^N 

+ iC' A [(iVtViV) • (VAT* x aiV) + (ViV+iV) • (JV^ x VAT) 

+ iC' h {N^N){VN^ ■ a x ViV) + iC' 6 (N*ffN) ■ (ViV f x ViV) 
+ (C^SikSji + C' 8 5uSkj + C' g 5ij8ki) 

+ (CVttfy + Cudfltfjy + C( 2 <%<y (iVta fc 9 i iV)(9 i iVta z iV) 



(C.l) 
(C.2) 
(C.3) 



(C.4) 
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Here we have shown explicitly only the operators TC K leading to nonvanishing r]H K rj, X 1 H K X 1 , 
\ l H K r], X 2 H K r], X 2 H K X 1 and h. c, which enter eqs. ( |4,9| )-( p~T3 ), The operators with three or more 



pion fields are irrelevant. The symbol means that the appropriate products in co-ordinate 
and isospin space have to be taken. 

Note that one has in principle four possible contact terms (two more involving t) in eq. (|C.3| ). 
However due to Fierz rearrangement they can be reduced to two. Further, the Hamiltonian used 
in this paper is always taken in normal ordering. 

The Hamiltonian given in contains apart from the terms enumerated above the two addi- 
tional interactions 



NUta- -Vtt ) V 2 N + V 2 iV f ( ta • Vtt J N 



+ -^V7V f f ta ■ -Vtt") • ViV 



(C.5) 



which lead to significant contributions to the two-nucleon potential at next-to-leading order. 
We now show that both terms (|C.5| ) can be completely eliminated from the Hamiltonian by 
integrating by parts and using the nucleons' equation of motion. For this let us consider the 
appropriate terms in the relativistic Lagrangian which contains exactly two nucleons, one pion 
and three derivatives. One can express them as 

(Na^P [bx*d v *E p + bfdfijtip + b 3 t^ v ~3 p + bit^tutp} tN) ■ 7T + h.c. . (C.6) 

Note that one can easily put all other terms with one or more derivatives acting on the pion 
fields into the form ( |C.G| ) by integrating by parts. The most general structure of the quantity 

a^" = ail Pg up + a 2 Yg^ + a 3 7 P <T + a^ l5 g^ + + (C.7) 

Now one can see, that all terms with the totally antisymmetric tensor e^ vpa do not contribute 
because of the derivatives commute (i.e. are symmetric under interchange of the indices) 

d p d v N = d v d^N . 

For all remaining terms with the metric tensor g^ v one can make use of the equation of motion 
of nucleon via 

(i 7 ^ -m + ...)N = (C.8) 

where the ellipsis represents terms of higher chiral dimension. We conclude that there are no 
corresponding terms in the relativistic Lagrangian. Nevertheless one could expect that terms 
of the type eq. ( |C.5| ) arise in the Lagrangian with fixed coupling constants after performing the 
nonrelativistic expansion for nucleons. In the one-nucleon sector such nonrelativistic expansion 
can be carried out explicitely using the path integral formulation heavy baryon formalism or 
by means of extended reparametrization invariance. One finds that one of the terms (|C.5| ) 
indeed enters the Lagrangian representing a 1/m 2 relativistic correction [Q]. However when 
the nucleon mass is counted via ( A.14| ), this term is obviously not relevant to the order we 



are working. Only the leading 1/m correction with exactly two derivatives resulting from the 
Lagrangian of lowest chiral dimension and representing the nucleon kinetic energy is to be taken 
into account explicitely. Our conclusion about the absence of terms ( |C.5| ) in the pion-nucleon 
Lagrangian agrees with corresponding result in the one-nucleon sector Hf^i . Finally it should 



be mentioned, that the structure ( C.5 ) can not arise by going from the Lagrangian to the 
corresponding Hamiltonian Q. 
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Figures 




1 2 3 

Figure 1: First corrections to the NN potential in the projection formalism: 
One-pion exchange and contact diagrams. The heavy dot denotes an insertion 
from the next-to-leading order pion-nucleon Lagrangian. Solid and dashed 
lines are nucleons and pions, respectively. 
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Figure 2: First corrections to the NN potential in the projection formalism: 
Two-pion exchange diagrams. Graphs which result from the interchange of 
the two nucleon lines are not shown. For notations, see fig. 1. 
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Figure 3: First corrections to the NN potential: Irreducible self- 
energy and vertex correction graphs. For notations, see fig. 1. 
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Figure 4: First corrections to the NN potential: Reducible self- 
energy graphs. For notations, see fig. 1. 
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Figure 5: Reducible graphs. In (a), a truly reducible diagram 
is shown. In the projection formalism, one has graphs like (b) 
and (c). These correspond to diagram 9 in fig. 2. The horizontal 
dashed lines count the free energies of the particles cut. 









Figure 6: One-loop corrections to the four-nucleon contact in- 
teraction. For notations, see fig. 1. 
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Figure 7: Leading diagrams contributing to the three-nucleon 
force that vanish. (1) is a representative for all graphs involving 
the 7T7riVA^-vertex. (2,3,4) all involve one four-fermion interac- 
tion. For notations, see fig. 1. 
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Figure 8: Leading diagrams contributing to the three-nucleon 
force. From the reducible graphs, only one representative is shown 
(as explained in fig. 5). For notations, see fig. 1. 
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